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ABSTRACT 

Aims. We present a detailed study of the radial distribution of the multiple populations identified in the Galactic globular cluster to 
Cen. 

Methods. We used both space-based images (ACS/WFC and WFPC2) and ground-based images (F0RS1@VLT and WFI@2.2m 
ESO telescopes) to map the cluster from the inner core to the outskirts (~20 arcmin). These data sets have been used to extract high- 
accuracy photometry for the construction of color-magnitude diagrams and astrometric positions of ~ 900 000 stars. 
Results. We find that in the inner ~2 core radii the blue main sequence (bMS) stars slightly dominate the red main sequence (rMS) 
in number At greater distances from the cluster center, the relative numbers of bMS stars with respect to rMS drop steeply, out to ~8 
arcmin, and then remain constant out to the limit of our observations. We also find that the dispersion of the Gaussian that best fits the 
color distribution within the bMS is significantly greater than the dispersion of the Gaussian that best fits the color distribution within 
the rMS. In addition, the relative number of intermediate-metallicity red-giant-branch stars which includes the progeny of the bMS) 
with respect to the metal-poor component (the progeny of the rMS) follows a trend similar to that of the main-sequence star-count 
ratio A'bMs/A'rMS- The most metal-rich component of the red-giant branch follows the same distribution as the intermediate-metallicity 
component. 

Conclusions. We briefly discuss the possible implications of the observed radial distribution of the different stellar components in 
(1) Cen. 

Key words. Globular clusters: general - Globular clusters: individual (w Cen [NGC 5139]) - Stars: evolution - Stars: Population 11 
-Techniques: photometric 

1 . Introduction w Cen hosts different stellar populations, most of them clearly 

visible in most of their evolutionary phases. 
^ The globular cluster (GC) Centauri is the most-studied stellar ^j^^^^ populations have been linked to the aforementioned 
H ; system of our Galaxy, but nevertheless one of the most puzzhng. j^et^nicity peaks, in photometric studies of the red-giant branch 
'^''^^ n'fnf^ "f''^^'''''^ * ^'a^f (RGB) (Pancino et al. |200i Hilker & Richder|200i Solhma 

im Norns & Bessell 1975 1977, R^E^" & dodgers 1975; EOOSit . the subgiant branch (SGB) (Hilker et al. |2004l 

^ell & Norris 1976^er et al. [M|Noms & Da Costa 3^^^^ |2005bl Stanford et al. « ViUanova et al. 

[19951 Suntzeff & KraftM Norris et al.JTMK with a primary 2007), and the main sequence (MS) (Piotto et al. 200?). The 
component at [Fe/H] ~ -1.7 to -1.8, and a long tail extending ^^^^ j.^^ ^^^^^^.^ ^ discovered by Piotto et 

up to [Fe/H] ~ -0 6 containing three or four secondary peaks ^^^^ ^^^^^^ ^^^^^ expectation from 

(see Johiison et al. 2009 for a recent update). It h^n shown, stellar-structure theory, the bluer of the two principal main se- 
both x^round-based photometry (Lee et aLpU Pancino ^^^^^^ ^^^^^ ^^^^ metal-rich than the redder one (rMS). 
et Rey et al ^ Sollima et al Mil Villanova V only possible way of reconciling the specti-oscopic observa- 

nt al. |2QQ3and Hubble Spacemescope (HST) P^metry ^.^^^ ^.^^ photometric ones is to assume a high overabun- 
(Anderson [19971 Bedin et al. [20041 FeiTaro et al. EOQl, that ^^^^^ ^j^^^ ^g^^.^ ^ Norris [20041 

Piotto et al. 12005 1 ). How such a high He content could have been 
formed is still a subject of debate (see Renzini 120081 for a re- 
Send offprint requests to: Bellini, A. view) 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which of the scenarios proposed to account for all the ob- 
is operated by AURA, Inc., under NASA contract NAS 5-26555, and served features of w Cen is a tidal stripping of an object that was 
on observations made with ESO telescopes at La Silla and Paranal originally much mor e mass ive (Zinnecker et al. 1988 ; Freeman 
Observatories. |1993i Dinescu et al. 119991 Ideta & Makino i2004t Tsuchiya et 
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al. (|2004] l; Bekki & Norris |2006l Villanova et al. |2007] ). In this 
scenario, the cluster was born as a dwarf elliptical galaxy, which 
was subsequently tidally disrupted by the Milky Way. Since all 
the populations of such a galaxy pass through the center, the nu- 
cleus would have been left with a mixture of all of them. 

It has also been suggested (Searle [T977l Makino et al. |1991| 
Ferraro et al. |2002| l that w Cen could have been formed by 
mergers of smaller stellar systems. In apparent support of this 
scenario, Ferraro et al. ( 120021) claimed that the most metal-rich 
RGB component of oj Cen (RGB-a, following the nomenclature 
of Pancino et al. I2000I I has a significantly different mean proper 
motion from that of the other RGB stars, and they concluded 
that RGB-a stars must have had an independent origin. However, 
Platais et al. ( 120031 ) showed that the proper-motion displacement 
seen could instead be an uncalibratable artifact of the plate solu- 
tion. More recently Bellini et al. (2009), with a new CCD-based 
proper-motion analysis, were able to demonstrate that all oj Cen 
RGB stars share the same mean motion to within a few km/sec. 
Anderson & van der Marel (2009 ) also find that the lower- turnoff 
population (the analog of the RGB-a) shows the same bulk mo- 
tion as the rest of the cluster Thus there is no longer a reason to 
think this population is kinematically distinct and an indication 
of a recent merger. Another indication that the cluster likely did 
not form by mergers can be found in the observation in Pancino 
et al. ( 120071 ) that all three RGB components share the cluster ro- 
tation, which would not be the case if different populations had 
different dynamical origins, or at least would require an unlikely 
degree of fine tuning. 

While OJ Cen was long thought to be the only cluster to ex- 
hibit a spread in abundances, we now know that it is not alone. 
M54 also clearly exhibits multiple RGBs (Sarajedini & Layden 
( IT9951 ); Siegel et al. 12007] ), SGBs (Piotto [2009] ), and has hints 
of multiple MSs. The complexity of M54 makes good sense, 
because it coincides with the nucleus of the tidally disrupting 
Sagittarius dwarf-spheroidal galaxy. M54 might be the actual 
nucleus or, more likely, it may represent a cluster that migrated 
to the nucleus as a result of dynamical friction (Bellazzini et 
al. 120081 ). a> Cen and M54 are the two most massive GCs in 
our Galaxy, and it is quite possible that they are the result of 
similar — and peculiar — evolutionary paths (Piotto [2009] ). In any 
case, even oj Cen and M54 are not the only clusters to exhibit 
non-singular populations. Exciting new discoveries, made in the 
last few years, clearly show that the GC multi-population zoo is 
quite populated, inhomogeneous, and complex. 

Piotto et al. (2007) published a color-magnitude diagram 
(CMD) of the globular cluster NGC 2808, in wh ich th ey iden- 
tified a well-defined triple MS (D'Antona et al. II2005II had al- 
ready suspected an anomalous broadening of the MS and had 
associated it with the three populations proposed by D'Antona 
& Caloi 1200411 to explain the complex horizontal branch (HB) 
of this cluster). Another globular cluster, NGC 1851, must have 
at least two distinct stellar populations. In this case, the obser- 
vational evidence comes from the split of the SGB (Milone et 
al. ,2008) . There are other GCs which undoubtedly show a split 
in the SGB, like NGC 6388 (Moretti et al. |2009] ) , M22 (Piotto 
12009] Marino et al. (i2009i) , 47 Tuc (Anderson et al. |2009] ), which 
also shows a MS broadening, or in the RGB, like M4 (Marino 
et al. 120081 ). Recent investigations (Rich et al. 2004'; Faria et al. 
120071 ) suggest that also other galaxies might host GCs with more 
than one population of stars. 

Multiple-population GCs offer observational evidence that 
challenges the traditional view. For half a century, a GC has 
been considered to be an assembly of stars that (quoting Renzini 
& Fusi Pecci II988I ): ''represent the purest and simplest stellar 



populations we can find in nature, as opposed to field popula- 
tions, which result from an admixture of ages and compositions'". 
If we allow for the fact that all the GCs for which Na and O 
abundances have been measured show a well defined Na/O anti- 
coiTelation (Carretta et al. 120061 120081) . suggesting an extended 
star-formation process, and that II of the 16 intermediate-age 
Large Magellanic Cloud GCs have been found to host multi- 
ple populations (Milone et al. |2009| ), multi-populations in GCs 
could be more the rule than the exception. De facto, a new era 
in globular-cluster research has started, and understanding how 
a multiple stellar system like o) Cen was born and has evolved is 
no longer the curious study of an anomaly, but rather may be a 
key to understanding basic star-formation processes. 

One way to understand how the multiple populations may 
have originated is to study the spatial distributions of the dif- 
ferent populations, which might retain information about where 
they formed. In particular, theoreticians have been finding that 
if the second generation of stars is formed from an interstellar 
medium polluted and shocked by the winds of the first genera- 
tion, then we would expect that the second generation would be 
more concentrated towards the center of the cluster than the first 
one ( see D' Ercole et al. 2008; Bekki & Mackey 2009; Decressin 
et al. 120081) . In the last of these references it is shown that in 
such a scenario the two generations of stars would interact dy- 
namically and would homogenize their radial distributions over 
time. As such, spatial gradients represent a fading fossil record 
of the cluster's dynamical history. 

Since oj Cen has such a long relaxation time (1.1 Gyr in the 
core and 10 Gyr at the half-mass radius, Harris [T996t , it is one of 
the few clusters where we might hope to infer the star-formation 
history by studying the internal kinematics and spatial distribu- 
tions of the constituent populations. These measurements will 
provide precious hints and constraints to allow theoreticians to 
develop more reliable GC dynamical models. 

In a recent paper, SoUima et al. (120071 ) showed that the star- 
count ratio NhMs/NiMS is flat beyond ~ 12', but that inward to 
~ 8' it increases to twice the envelope value. Thus the bMS stars 
(i.e., the supposed "He-enriched" population) are more concen- 
trated towards the center than the rMS, which is presumed to be 
the first generation. Unfortunately, Sollima et al. ( 120071 could 
not provide information about the trend of A^bMS /A^rMS within 
~ 8', which corresponds roughly to 2 half-mass radii (r/,). 

On the other hand, the radial distribution of RGB subpopu- 
lations has been analyzed by many authors (Norris et al. 1997] 
Hilker & Richtler 12^0] P ancino et al. 2000 2003; Rey et al. 
2004i SolUma et al. l2005al CasteUani et al. 2007; Johnson et al. 
120091 ). All these works agree that the intermediate-metallicity 
population (RGB-Mint) is more centrally concentrated than the 
more metal-poor one (RGB-MP). However, there is a disagree- 
ment about the most metal-rich population (RGB-a): Pancino et 
al. (2000), Norris et al. (1997), and Johnson et al. (2009) found 
that the most metal-rich stars (RGB-a) are as concentrated as the 
intermediate-metallicity ones, and consequently more concen- 
trated than the most metal-poor stars, whereas Hilker & Richtler 
(120001 ) and Castellani et al. (120071 ) considered the RGB-a com- 
ponent to be the least-concentrated population. (Since our work 
in progress was already favoring the former view over the latter, 
we were concerned to reach the definitive truth of this matter). 

In the present paper, we trace the radial distribution of the 
stars of OJ Cen , both on the MS and in the RGB region. Our ra- 
dial density analysis covers both the center and the outskirts of 
the cluster, taking advantage of the combination of four instru- 
ments on three different telescopes, and of our proper-motion 
measurements on ground-based multi-epoch wide-field images 
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(Bellini et al. |2009l l. In Section 2 we describe in detail the photo- 
metric data and the reduction procedures. Section 3 presents our 
analysis of the radial distribution of the stars on the two MSs. In 
Section 4 we perform an analogous study for the RGB stars. A 
brief discussion follows in Section 5. 

2. Observations and data reductions 

To trace the radial distribution of the different stellar populations 
in a> Cen, we analyzed several data sets, from four different cam- 
eras. To probe the dense inner regions of the cluster we took ad- 
vantage of the space-based high resolving power of HST, using 
both the Wide Field Channel (WFC) of the Advanced Camera 
for Surveys (ACS), and the Wide Field and Planetary Camera 2 
(WFPC2). For the relatively sparse outskirts of the cluster, we in- 
stead made use of deep archival ground-based observations col- 
lected with the FORS 1 camera of the ESO Very Large Telescope 
(VLT). In addition, to link all the different data sets into a com- 
mon astrometric and photometric reference system, we used the 
Wide Field Imager (WFl) at the focus of the ESO 2.2m telescope 
(hereafter WFl@2.2m). This shallower data set was also used to 
study the red-giant branch in the outskirts of the cluster 

Figure [T] shows the footprints of the data sets, centered 
on the recently determined accurate center of a> Cen: RA - 
13:26:47.2 4, Dec = -47:28:46.45 (J2000.0, Anderson & van 
der Marel I2009I I. The red footprints are those of HST obser- 
vations. The larger ones are the ACS/WFC data sets, a 3 x 3 
mosaic centered on the cluster center and a single field ~17' 
SW of the center. The smaller red field, ~7' S of the center, 
was observed with WFPC2. Blue rectangles show the partially 
overlapping FORSl@VLT fields, extending from ~6' to ~25'. 
The large field in magenta is the ~33' x 33' field-of-view of our 
WFl@2.2m proper-motion catalog (Bellini et al. 120091 ). The fig- 
ure also shows the major and minor axes (solid lines), taken from 
van de Ven et al. ( 2006 ). We divided the field into four quadrants, 
centered on the major and minor axes. The quadrants are labeled 
with Roman numerals and separated by dashed lines. We will use 
them to derive internal estimates of the errors of the star-count 
distribution. Concentric ellipses, aligned with the major/minor 
axes, have ellipticity of 0.17, coincident with the average ellip- 
ticity of (L> Cen Geyer et al. 119831 . These ellipses will be used to 
define radial annuli, in Section 2.8. Thick black circles mark the 
core radius (re - I'A) and the half-mass radius (r/, = 4fl8) (from 
Harris [1996). If we assuming a cluster distance of 4.7 kpc (van 
de Ven et al. 2006; van der Marel & Anderson |2009l l, the two 
radii correspond to 1.9 pc and 5.7 pc, respectively. 

The details of the data sets are summarized in Table [T] In the 
following subsections we give brief descriptions of the reduction 
procedures, which have been presented in more detail in various 
other papers. The FORSl data, however, were taken by Sollima 
et al. ( I2007I I, for a purpose similar to ours; we will give a full 
description of our reduction in subsection |2.4| 

2. 1 . HST: ACS/WFC inner 3x3 mosaic 

This data set (inner nine red rectangles in Fig. [T] GO-9442, PI 
A. Cool) consists of a mosaic of 3 x 3 fields obtained with the 
ACS/WFC through the F435W and F625W filters. This camera 
has a pixel size of ~50 mas and a field of view of 3f3 x 3'3. Each 
of these nine fields has one short and three long exposures in 
both F435W and F625W. The mosaic covers the inner -lO'xlO', 
the most crowded region of co Cen. These images, which were 
used by Ferraro et al. (2004) and by Freyhammer et al. (2005), 
and which we used in both Bedin et al. (,2004), and Villanova et 
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Fig. 1. The footprints of the ll) Cen data sets used in this work. 
North is up, east to the left. The "*" marks the cluster center The 
3x3 ACS/WFC mosaic (in red) is that of GO-9442, while the 8 
FORSl fields are marked in blue. The largest field (in magenta) 
comes from the WFl proper-motion catalog (Bellini et al. 2009). 
This wide-field catalog has been used to register the FORS 1 and 
the ACS/WFC inner-mosaic data into a common astrometric and 
photometric reference system (see text). The smaller WFPC2 
field at ~T south and the outer ACS/WFC field at -17' from 
the cluster center are also shown (in red). Concentric ellipses, 
centered on the center of ll) Cen and aligned with the major and 
minor axes, show the radial bins that we created. Ellipses are 
split into quadrants (dashed lines), labeled with Roman numer- 
als. Each quadrant is bisected by the major or minor axis. Thick 
black circles mark the core radius (re - 1 '4) and the half-mass 
radius (r,, = 4:i8) (from Harris [T99^ . 



al. (2007), were reduced using img2xym_WFC.09x1 0, which is 
a publicly available FORTRAN program, described in Anderson 
& King ( I2OO6I 1. The program finds and measures each star 
in each exposure by fitting a spatially-variable effective point- 
spread function. The independent measurements of the stars 
were collated into a master star list that covers the entire 3x3 
mosaic field. For each star we constructed an average magnitude 
in each band, and computed the rms deviation of the multiple 
measurements about this average. Instrumental magnitudes were 
transformed into the ACS Vega-mag flight system following the 
procedure given in Bedin et al. ( 12005) , using the zero points of 
Sirianni et al. (I2005I I. Since the zero points are valid only for 
fluxes in the _drz exposures, we computed calibrated photome- 
try for a few isolated stars in the _drz exposures and used this 
to set the zero points for the photometry that was based on the 
individual _f It images. Saturated stars in short exposures were 
treated as described in Section 8.1 in Anderson et al. (2008)1. 
Collecting photoelectrons along the bleeding columns allowed 
us to measure magnitudes of saturated stars up to 3.5 magni- 
tudes above saturation (i.e., up to mF435w-12 mag), with errors 
of only a few percent (Gilliland 2004). We used the final cata- 
log, which contains more than 760 000 stars, to trace the radial 
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Table 1. Data sets used in this work. 



Data set 


Filter 


# images x Exp. time (s) 


3x3 ACS/WFC 


F435W 


27 X 340, 9x 12 


inner mosaic 








F625W 


27 X 340, 9x8 


ACS/WFC 


F606W 


2 X 1285, 2 x 1300, 


-17' 




2x 1331, 2 X 1375 




F814W 


4 X 1331, 2 X 1340, 2 X 1375 


WFPC2@ HST 


F606W 


2 X 300, 1 X 600 


~T 


F814W 


2 X 400, 1 X 1000 


FORSl@VLT 


B 


20 X 1100 




R 


20 X 395 


WFI@2.2m 


B 


1 X 10, 1 X 15, 11 x30, 1 x40. 






1 X 60, 1 X 120,2 X 240, 2 x 300 




Rc 


1 X 5, 1 X 10, 1 X 15, 1 X 30, 5 X 60 




V 


6 X 5, 9 X 10, 1 X 15, 3 X 20, 






2 X 30, 10 X 40, 4 X 45, 3 x 60, 






10 X 90, 7 X 120, 1 X 150, 3 x 240 



distribution of RGB and MS stars in this most crowded region 
of the cluster. 

2.2. HST: ACS/WFC outer field 

The outer ACS field (-17' SW of the cluster center, see Fig.[Tll 
comes from proposals GO-9444 and GO- 10101 (both with PI I. 
R. King), using the F606W and F814W filters. The photometry 
from the first-epoch observations was published in Bedin et al. 
( I2004I I. The photometry presented in the present paper comes 
from the full two-epoch data set for this field; the two epochs 
also allow us to derive proper motions and perform a critical 
cluster/field separation. A detailed description of the data reduc- 
tion, the proper-motion measurement, and the resulting CMDs 
will be presented in a forthcoming paper. The reduction and cal- 
ibration of these data sets use procedures similar to those used 
for the central mosaic, and provided photometry for ~3500 stars. 

2.3. HST: WFPC2 field 

We also make use of one WFPC2 field, ~1' south of the cluster 
center (see Fig.[T]). This data set consists of 2 x 300 4- 600s expo- 
sures in F606W, and 2 x 400 -i- 1000s in F8 14W (GO-5370, PI R. 
Griffiths), and contains 9214 stars. These images have been re- 
duced with the algorithms described in Anderson & King (2000). 
The field was calibrated to the photometric Vega-mag flight sys- 
tem of WFPC2 according to the prescriptions of Holtzman et 
al. ( 119951 ). This WFPC2 field is particularly important in tracing 
the distribution of stars in the MS of cj Cen, because it is at a 
radial distance from the center of the cluster where there are no 
suitable ACS/WFC observations and where ground-based obser- 
vations are almost useless because of crowding. 

2.4. VLT: eight F0RS1 fields 

The VLT data set consists of eight partially overlapping FORS 1 
fields, each with a pixel size of 200 mas and a field of view of 



multiple stellar populations in oj Centauri 



18 20 22 24 26 




IB 20 22 24 



B 



16 18 20 22 24 26 

Fig. 2. Selection criteria used to isolate FORS 1@ VLT stars for 
our MS subpopulation analysis. Panel (a) shows sharp values 
versus B magnitude, and panel (b) x versus B. Panels (c) and (d) 
show the photometric errors as a function of the B and Cousins- 
Rc magnitudes respectively. Only stars that passed the sharp 
selection criterion (black in the first panel), are plotted in the 
subsequent panels; similarly, only stars that also survived the;^f 
selection are shown in the remaining two panels. 

6^8 x 6^8. These fields (the blue rectangles in Fig.[T]) probe the 
regions between 6' and 25' from the center of a> Cen . The set of 
images consists of 20 x 1 100s exposures in B, and 20 x 395s in 
R, and are the same images used by SoUima et al. ("2007). There 
are four images in each field (two per filter), except that the third 
and fourth fields have four images per filter (see Fig. [1] for field 
numbers). This is the only data set that we reduced specifically 
for the present work. For this reason we give a more detailed 
description of our reduction procedure. 

We retrieved the data sets from the ESO archive; master- 
bias and flat-field frames were constructed using standard IRAF 
routines. Photometric reduction of the images was performed 
using P Stet son's DAOPHOT-ALLSTAR-ALLFRAME pack- 
ages (Stetson ' ^198711 199 4). For each exposure we constructed a 
quadratic spatially variable point-spread function (PSF) by using 
a Penny functiorQ, and for each individual exposure we chose — 
by visual inspection — the best 100 (at least) isolated, bright, un- 

' A Penny function is the sum of a Gaussian and a Lorentz function. 
In this case we used five free parameters: half- width at half-maximum 
of the Penny function, in the x and in the y coordinate; the fractional 
amplitude of the Gaussian function at the peak of the stellar profile; the 
position angle of the tilted elliptical Gaussian; and a tilt of the Lorentz 
function in a different direction from the Gaussian. The Lorentz func- 
tion may be elongated too, but its long axis is parallel to the x oi y 
direction. 
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saturated stars that were suitable for mapping the PSF variations 
all over the image. We used ALLFRAME on each individual 
field, keeping only stars measured in at least four images. The 
photometric zero points of each field were registered to the in- 
strumental magnitudes of the fourth field (the less crowded of 
the two that have more exposures). Finally, photometric and 
astrometric calibration was performed using the WFI@2.2m 
astro metric-photometric catalog by Bellini et al. (2009) as a ref- 
erence. As a result, we brought the FORSl R magnitudes to the 
Cousins-7?c photometric system used by WFI@2.2m. Our final 
FORSl catalog contains -133 000 objects. 

Since the innermost FORSl field is seriously affected by 
crowding, we did not use it in the present analysis. Fig.|2]plots 
the sharp, x, and as and ctr^ calculated by ALLFRAME, as 
functions of stellar magnitude, for the stars in the FORS 1 cat- 
alog. To choose the well-measured stars, we drew by eye the 
cut-off boundaries in the quality parameters that retained objects 
that were most likely to be well-measured stars. Panel (a) shows 
sharp values versus B magnitude. Stars that passed the selection 
criterion are shown in black. Panel (b), which includes only stars 
that passed the sharp cut, shows x values versus B. Stars that 
also passed the x criterion are in black. In panel (c) we plot the 
(Tb values versus B, for the stars that survived these two selec- 
tions. Again, the stars with good photometry are shown in black. 
Finally, in the last panel we plot cr^^ values versus Rc, for all the 
survivors, and we highlight in black those that survived this se- 
lection too. At the end of these selection procedures, we are left 
with a catalog of ~66 500 stars. We note that while these selec- 
tion criteria affect stars at different magnitudes differently, they 
should not affect the ratio of stars on the bMS and rMS, since at 
a given magnitude the two populations should both have about 
the same photometric error, and the same probability of making 
it into our catalog. 

2.5. WFI@2.2m 

This data set was collected at the 2.2m ESO Telescope, with 
the WFI camera, between 1999 and 2003. The WFI@2.2m cam- 
era is made up of a mosaic of 4 x 2 chips, 2048 x 4096 pixels 
each, with a pixel scale of 238 mas/pixel). Thus, each WFI expo- 
sure covers ~34' x 33'. The ai Cen astrometric, photometric, and 
proper-motion catalog based on this data set and presented in 
Bellini et al. (|2009| is public, and contains several wide-band 
(U,B,V,Rc,Ic) filters plus a narrow-band filter (658nm), and 
covers an area of ~33' x 33' centered on the cluster center. We 
refer the reader to Bellini et al. (2009^ for a detailed discussion 
of the data-reduction and calibration procedures. 

Briefly, photometry and astrometry were extracted with the 
procedures and codes described in Anderson et al. ( 120061 ) . 
Photometric measurements were corrected for "sky concentra- 
tion" effect^ and for differential reddening, as described in 
Manfroid & Selman (2001 ) and BelHni et al. ( |2009l l. Global star 
positions are measured to better than ~45 mas in each coordi- 
nate. Photometric calibration in the B, V, Rc, Ic bands is based 
on a set of ~3000 secondary standard stars in to Cen , avail- 
able on-line (Stetson 200(T, 2005 ). Color equations were derived 
to transform our instrumental photometry into the photometri- 
cally calibrated system using an iterative least-squares linear fit. 
Thanks to the four-year time-baseline, we were able to success- 
fully separate cluster members from field stars by means of the 
local-transformation approach (Anderson et al. i20061 l. giving us 



^ Light contamination caused by internal reflections of light in the 
optics, causing a redistribution of light in the focal plane. 




0.5 1 1.5 a 



B-Rc 

Fig. 3. {Top left:) Selected stars in common between the 
ACS/WFC 3x3 mosaic and the WFI@2.2m data sets. {Top 
right:) Horizontal lines mark the magnitude interval used to de- 
rive calibration equations. {Bottom panels:) Calibration fits used 
to transform Vega-mag ACS/WFC mF435w and mpeasw magni- 
tudes into the WFI@2.2m photometric system. See text for de- 
tails. 



proper motions more precise than ~ 4 mas yr down to B ~20 
mag, for ~54000 stars. 

2.6. The astrometric and photometric reference frame 

The large field of view of the WFI@2.2m camera makes our 
WFI catalog an ideal photometric and astrometric reference 
frame to which to refer all the other observations, from different 
telescope-camera-filter combinations. For each catalog we made 
the tie-in by means of stars that were in common with the refer- 
ence catalog. For positions we derived a general six-parameter 
linear transformation to the astrometric system of the WFI cat- 
alog. For photometry we used as a reference standard the B and 
Cousins-7?c magnitudes of the WFI@2.2m catalog, and trans- 
formed the magnitudes of each other catalog to this standard. For 
the mF435w and mF625w magnitudes of the central mosaic of 3x3 
ACS/WFC fields, we used ~3300 stars that had been observed 
in common, located outside 4' from the cluster center to avoid 
the most crowded regions in the WFI data set (top-left panel of 
Fig.|3]). We excluded from this sample saturated stars in the WFI 
data set, keeping only the brighter (14.9 < B < 16.5) and well 
measured {<Jb,Rc < 0.02 mag) ones (top-right panel in Fig. [3). 
The adopted calibration fits are shown in the bottom panels of 
Fig.|3] We did similarly for the FORSl B and R magnitudes. 
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Fig. 4. B vs. B - Rc CMD of WFI@2.2ot stars, after calibration 
and proper-motion selection (see Bellini et al. ,2009j . 



ACS/WFC mosaic. It is clear that the MS population can be stud- 
ied in all but the WFI CMD, and the RGB population can be 
studied in the WFI and inner ACS data sets. 

Now that we have a comprehensive sample of oj Cen stars, 
both for the bright stars and for the faint ones, covering the cen- 
tral region all the way out to ~25', we can define robust selection 
criteria for the subpopulations to track how the population frac- 
tions vary with radius. 

2.8. The angular radial distance: r* 

Since w Cen is elongated in the plane of the sky, it does not 
make sense to analyze its radial profile via circular annuli. We 
therefore decided to extract radial bins in the following way. 
We adopted the position angle (P. A.) of 100° for the major- 
axis (van de Ven et al. .2006 ), and an average ellipticity of 0.17 
Geyer et al. II983I I. To define the bins of the radial distribution 
we adopted elliptical annuli, whose major axes are aligned with 
the o) Cen major axis, and stars were extracted accordingly (see 
Fig. [TJ. To indicate the angular radial distance from the cluster 
center, we used the equivalent radius r*, defined as the radius of 
the circle with the same area as the corresponding ellipse (i.e., 
the geometrical mean of the semi-major and semi-minor axes). 
Each of the small fields (the outer ACS field and the WFPC2 
field), we considered as a single radial bin. 



Calamida et al. (I2005I I measured a diff'erential reddening of 
uptoE(B-y) ~0.14in aregionof ~14'x 14' centered on w Cen. 
This result has been questioned by Villanova et al. (120071 ); in 
their Figs. 1-6, the sharpness of the SGB sequences suggests 
that the existence of any serious differential reddening is very 
unlikely. But in any case, a proper radial-distribution analysis 
needs coiTection even for a differential reddening that is of the 
order of few hundredths of a magnitude. Our coiTections for dif- 
ferential reddening followed the method outlined by Sarajedini 
et al. ( I2007I I. which uses the displacements of individual stars 
from a fiducial sequence to derive a reddening map. 

The outer ACS/WFC field at ~17' from the cluster center 
and the WFPC2 field at ~7' provide stellar photometry in the 
F606W and F814W bands. For the ACS field we have overlap 
with the WFI catalog, which allows us to calibrate the photom- 
etry, but the stars available are all on the main sequence above 
'WF606W - 21, so they have a very narrow range in color, and we 
cannot empirically determine the color term in the calibration. 
For the WFPC2 field, in addition to the problem of the limited 
color baseline, the WFI photometry in this inner field is of low 
quality on account of ground-based crowding. For these reasons, 
we decided to not transform the photometry of these two fields 
into the photometric reference system of WFI@2.2m, but dealt 
with them in the HST Vega-mag flight system. 

2. 7. The deep color-magnitude diagrams 

Our proper-motion-selected WFI@2.2m B \s. B - Rc CMD is 
shown in Fig.|4] All the main features of the cluster are clearly 
visible, except for the split MS, since the WFI data go down only 
a magnitude or so below the turnoff. The CMDs of the other 
data sets that we analyzed are presented in Fig.|5] where the top- 
left panel refers to the eight FORSl@VLT fields, the middle- 
left panel to the proper-motion-selected CMD of the external 
ACS/WFC, the bottom-left panel the CMD from the WFPC2 
field, and the right panel of Fig. |5] the CMD of the inner 3x3 



3. MS subpopulations 

Our goal in putting together these varied catalogs is to quantify 
the differences in the radial profiles of the various subpopula- 
tions of CO Cen. One way to do this would be to measure the sur- 
face density profile for each group and compare them directly, 
but this would require accurate completeness corrections and 
careful attention to magnitude bins. Since our interest, however, 
is simply to determine how the populations vary with respect to 
each other, we need only measure the ratio of the populations as 
a function of radius. This ratio should be independent of com- 
pleteness coiTections and the details of the magnitude bins used, 
since the bMS and rMS differ only slightly in color and are ob- 
served over the same magnitude range. 

Our analysis of the NhMs/Nrus ratio is based mostly on 
the data sets from the inner ACS/WFC 3x3 mosaic and 
FORSl@VLT, which allow us to map the ratio of bMS/rMS 
from the cluster center out to ~25', once the photometry and 
astrometry have been brought into the same reference system. 
The other two fields, each of which covers only a small region, 
provide only one point each in our analysis of A^bMS /Nms versus 
radius. Moreover, since we were not able to bring ;7Jf606w and 
'WF814W photometry of the outer ACS and the WFPC2 field into 
the WFI B and Rc photometric system, we kept the WFPC2 and 
the outer ACS/WFC data sets in their native photometric sys- 
tem, and used them only for a further (though important) con- 
firmation of the radial gradient found with the FORS 1 and inner 
ACS/WFC data sets. 

3.1. Straightened main sequences 

In order to analyze the color distribution of the stars along the 
MS in a more convenient coordinate system, we adopted a tech- 
nique previously used with success in w Cen (Anderson 1 19971 
2002), and in other works (Sollima et al. 120071 Villanova et al. 
2007; Piotto et al. l2007l Milone et al. l2008l Anderson et al. 
,2009i l. 
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ACS/WFC 3x3 




Fig. 5. (7b;? left): CMD from the eight FORSl @VLT fields. We can measure stars from the bottom of the RGB down to B -27.5 
mag. {Middle left): proper-motion-selected CMD from the outer ACS/WFC field. {Bottom left): CMD from the WFPC2 images 
located ~7' south of the cluster center. {Right panel): CMD of the inner 3x3 ACS/WFC fields. In the top left and the right-hand 
CMDs, the bMS and rMS fail to show separately only because the profusion of points blackens their whole region. 



We defined fiducial fines in the CMDs (drawn by hand), such 
as to be equidistant from the ridge lines of the bMS and rMS 
stars. We avoided choosing the ridge line of either sequence as 
our fiducial fine, because we wanted a system in which both 



the sequences are as parallel and as rectified as possible. We 
used different fiducial lines for the B, B - Rq CMDs of the in- 
ner ACS/WFC and the FORSl data sets and for the (otf606w, 
OTF606W - OTF814W) CMDs of the WFPC2 and outer ACS/WFC 
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Fig. 6. The left panel shows a randomly selected 8% of the stars 
in the CMD of the inner -10' x 10' ACS/WFC images, in the 
region of the MS where the two branches are most separated in 
color. The middle panel shows the CMD of the F0RS1@VLT 
fields. The right panels show the outer ACS/WFC field (bottom) 
and WFPC2 field (top). The MS duality is clearly detected in 
all diagrams (see also Fig.|7]i. The dashed horizontal lines mark 
the selected magnitude range for the definition of the bMS and 
rMS samples used in the derivation of their radial profiles. The 
fiducial lines (drawn by hand) that were used to straighten and 
separate the sequences are also plotted (in red in the color ver- 
sion). 



data sets. In this way, we were sure to straighten the MSs in the 
same consistent way for the two different sets of filters. Then we 
subtracted from the color of each star the color of the fiducial 
line at the same luminosity as the star 

In Fig.|6]we show the CMDs in the oj Cen MS region for the 
central mosaic of ACSAVFC data (left panel), the FORSl @VLT 
(middle panel), and the WFPC2 ~T field and the ACS/WFC 
field at -IT (right panels). In the case of the central ACS/WFC 
data, we plotted only a randomly chosen 8% of the stars, in 
order to show the two sequences clearly. In all the CMDs the 
MS splitting is clearly visible. For the inner ACS/WFC and 
FORS 1 data sets we restricted our MS analysis to the magnitude 
range 20.9 < B < 22.1 (dashed lines in Fig.|6]l, the interval in 
which the two MSs are most separated in color and are parallel. 
For the same reasons we analyzed stars in the magnitude range 
20.6 < mF606w < 21.9 for the WFPC2 and the outer ACS/WFC 
data sets. The bright limit also avoids the saturated stars in the 
deep WFC exposures. The adopted fiducial lines are again plot- 
ted (in red in the color version of the paper). 

In Fig.|7]we show straightened CMDs for the same data sets 
shown in Fig.|6l with the only difference being that we now plot a 
20% randomly generated sample of stars for the inner ACS/WFC 
data set, since the expanded color baseline allows more points to 
be seen. It is worth noting that even a simple inspection shows 
the A^bMS /NrM5 ratio clearly decreasing as we go from the central 
cluster regions to the outer ones. It is also clear that the spread 
in the bMS is somewhat greater than that of the rMS. 



ACS/WFC inner 
3x3 fields 
20% randomly 
selected stars. . 



J.l.' 
0.05 



-Acs/wi;c~i:j'--. ■ ■ ■ . - 

0.1 -0.1-0.05 0.05 0.1 -0.1-0.05 0.05 0.1 



A{B-R,) Aim,, 
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Fig. 7. Same as Fig.|6] but after subtraction, from the color of 
each star, of the color of the fiducial line at the same luminosity. 
In the left panel we show a randomly selected 20% of the stars 
from the ACS/WFC central-mosaic data (rather than the previ- 
ous 8%, since the color-scale is now less compressed). 



Finally, note that we call the color deviation of a star from 
the fiducial line A(B - Re). We shall use this notation frequently 
in what follows. 

Our aim in selecting the best-measured stars in the previ- 
ous sections was so that we would be able to assign the stars to 
the different populations as accurately as possible. Similarly, as 
much as possible we transformed our photometry into the same 
system, so that our population selections throughout the cluster 
would be as consistent as possible. 

Even with these careful steps, however, it is still difficult to 
ensure that we are selecting stars of the same population in the 
inner parts of the cluster as in the outer parts. Even if we had ob- 
servations with the same detector at all radii, the greater crowd- 
ing at the center would increase the errors there. On the other 
hand, our use of ground-based images for the outer fields actu- 
ally makes those fields even more vulnerable to crowding effects. 

Another complication comes from main-sequence binaries, 
which at the distance of a globular cluster are unresolved. 
Relaxation, causing mass segregation, will concentrate them to 
the cluster center and cause a redward distortion of the main se- 
quence there. 

Moreover, in the lower-density outer regions of the cluster 
we can get the same statistical significance only by using larger 
areas, with an increased vulnerability to inclusion of field stars. 
Finally, the red side of the main sequence is contaminated by 
the anomalous metal-rich population (hereafter MS-a), which is 
clearly connected with RGB-a. Even if these stars include only 
~5% of the total cluster members (Lee et al. 1999^; Pancino et 
al. |2000 l SolUma et al. 2005a, Villanova et al. 2007), they are an 
additional source of pollution for rMS stars — against which we 
now take specific precautions. 
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Fig. 8. (Left panel): B versus A(B-Rc) diagram for selected stars 
in our data set from the ACS/WFC 3x3 mosaic, divided into five 
magnitude intervals. We now plot all the stars in this sample, 
not just a randomly selected subset. (Right panels): A(B - Rc) 
histograms with the Gaussian best fits. See text for details. 



3.2. Dual-Gaussian fitting 

There is no way of dealing with the above issues perfectly, but 
we did our best to make our measurements as insensitive to 
them as possible. To this end, we measured the bMS and rMS 
fractions by simultaneously fitting the straightened color dis- 
tributions with two Gaussians, and taking the area under each 
Gaussian as our estimate of the number of stars in each pop- 
ulation. By keeping the width of each Gaussian an adjustable 
parameter, we allowed in a natural way for the fact that the pho- 
tometric scatter differs from one radius and data set to another 

While the dual Gaussians provide a natural way of measuring 
the two populations in data sets that have different color base- 
lines and different photometric errors, there is one serious com- 
plication. As we have indicated, there is an unresolved, broad 
population of stars redward of the rMS that consists of blends, 
binaries, and members of the MS-a branch. Since it is unclear 
what relation this mixed population has with the two populations 
that we are studying, we wanted to exclude it from the analysis 
as much as possible. We did so by cutting off the reddest part of 
the color range, and confining our fitting to the color range that 
is least disturbed by the contaminated red tail. 

In order to choose the red cutoff as well as possible, we gath- 
ered together all of the stars in each data set. Below we will de- 
scribe for simplicity only the case of the central 3x3 mosaic of 
ACS images in B and Rc. The procedure followed is, however, 
the same for the other data sets. 

Within this data set we chose the MS stars that were in the 
magnitude range 20.9 < B < 22.1 (within which the two MSs 
are almost parallel and are maximally separated in color) and in 
the color range -0.25 < A(B - Rc) < 0.25 mag. We emphasize 
that this ensemble of the data set, within which we will later see 
a considerable gradient in the relative numbers of bMS and rMS 
stars, will not be used to derive population results in the case of 
the inner ACS/WFC data set, but only to choose the red cutoff. 



We divided these stars into five magnitude intervals, because the 
observational errors, which increase the spread of the sequences, 
depend on magnitude. Next, we plotted histograms of the A(B - 
Rc) distribution within each magnitude interval, using a bin size 
of 0.006 mag. This size is -1/4 of the typical photometric error 
in color; it makes a good compromise between a fine enough 
color resolution, on the one hand, and adequate statistics, on the 
other hand. 

The actual choice of the red cutoff is a two-tiered proce- 
dure. We must first develop a procedure for the fitting of dual 
Gaussians to a set of bins that has a red cutoff; then we must 
decide on a value of Med, the number of bins that we include on 
the redder side of the red Gaussian. 

Although from a mere inspection of the histograms it is clear 
where, approximately, the peak of the red Gaussian should lie, 
the narrowness of the bins leaves it uncertain in which particular 
bin the peak of the red Gaussian will actually fall. Since the red 
cutoff, Med, is defined as being counted from that bin, we had 
to resort to an iterative procedure to locate the cutoff for a given 
value of Med. We began by choosing a cutoff safely to the red 
of where we guessed that the cutoff would actually fall, and then 
using that cutoff in a first try at fitting the dual Gaussians. The 
iteration then consisted of placing the cutoff just beyond Med 
bins on the red side of the peak of the red Gaussian and fitting 
again; this new fit might cause the red peak to move to a different 
bin. When the red peak stays in the same bin, the iteration has 
converged; this happened after very few iterations. 

We assumed trial values of Mred from 2 to 5, and for each of 
those values we iteratively computed the Gaussian parameters 
for each of the five magnitude intervals. We chose as the best 
value for Mred the one for which the five values of MbMS /Mi-ms 
were the most consistent. This value turned out to be Mred - 3. 
With this choice made, we then moved on to fit dual Gaussians 
to each of our detailed data sets. 

Fig.[8]shows the results of this procedure. In the left panel we 
show our selected stars in the B versus A(B - Rc) diagram — all 
of the stars this time, rather than a random selection of a frac- 
tion of them. The horizontal lines delineate our five magnitude 
intervals. On the right we show the final A(B-Rc) histogram for 
each magnitude interval, and the dual-Gaussian fit to it. The in- 
dividual Gaussians are shown in blue and red, respectively, and 
the black curve is their sum. The vertical blue and red lines are 
the centers of the respective Gaussians, and the vertical black 
line shows the red cutoff. Note that we do not show the verti- 
cal boundaries between the bins of a histogram, because on this 
scale they would be too close to each other. Nor do we show 
the Poisson errors of the counts in the bins, because they are 
small and would obscure the bin values themselves; the size of 
the errors is amply clear from the smoothness of the values in 
neighboring bins. The counts in the histograms are normalized 
so as to make the height of the red Gaussian equal to unity. 

3.3. Tiie Radial Gradient of MbMs/M-MS 

Having chosen the position of the red cutoff, we were able to 
perform dual-Gaussian fitting on each of our data sets. Figure |9] 
shows our fits. We divided the inner ACS/WFC mosaic and the 
outer FORSl @VLT data sets into five radial intervals for each. 
The intervals were chosen in such a way as to have the same 
number of selected stars in each of them, so that the statistical 
sampling errors will be uniform. (The reader should note that 
Fig.|9]shows all of our fields, in radial order, so that the WFPC2 
field follows the inner ACS fields, and the outermost ACS field 
falls between two of the FORSl fields.) 
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Fig. 9. Dual-Gaussians fits. As in Fig. [8] the Gaussian fits to the bMS and rMS are in blue and red respectively, and their sum in 
black. The vertical dashed lines mark the centers of the individual Gaussians. The individual panels are arranged in order of effective 
radius. (Note that all our fields are shown here, in radial order, so that the WFPC2 field follows the inner ACS fields, and the outer 
ACS field falls between two of the FORSl fields.) 
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Fig. 10. A^bMS /NiMS ratio versus equivalent radius r*. Different colors and symbols refer to different data sets. Dashed vertical lines 
mark the core radius and the half-mass radius. Error bars were calculated from the residuals of values in individual subdivisions 
(quadrants for the inner ACS/WFC mosaic, magnitude intervals in each outer field). To improve the radial resolution for the out- 
ermost annulus of the inner ACS/WFC mosaic, we also divided it into four sub-annuli (crossed open circles). See text for a fuller 
explanation. 



Figure [10] shows our results for the radial variation of the field, and the outer ACS field. Symbols of a different shape dis- 
bMS to rMS ratio, for the five radial parts of the inner ACS mo- tinguish the various types of field. The outermost radial interval 
saic, the five radial intervals of our FORS 1 fields, the WFPC2 of the ACS/WFC mosaic is a special case, however, since it con- 
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Table 2. Dual-Gaussian fitting results. For each data set (first column) we give in Cols. 2-4 the radial extent (minimum, median, 
and maximum r*). In Cols. 5-8 are the sigmas of the Gaussians that best fit the bMS and rMS color distributions, with errors. In the 
next two columns are the A^bMS /A^rMS ratio and its error. The next-to-last column gives the color difference between the two Gaussian 
peaks, and the final column identifies the color baseline of the data set. 
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sists largely of the four corners of the mosaic, and it spans a 
larger radial extension. To better map the bMS/rMS distribution 
in this radial interval, we decided to further split it into four sub- 
annuli. In this way we increase the radial resolution, but pay the 
price of larger sampling errors. We have therefore plotted the 
outermost radial interval of the inner ACS/WFC mosaic twice, 
once as a whole annulus, and once as four sub-annuli (marked 
as crossed open circles in Fig.fTOll. 

Our choice of using ellipses with fixed ellipticity and posi- 
tion angle to extract radial bins could have introduced some sys- 
tematics in our derived A^bMs/M-MS ratios. To address this issue, 
we recalculated the A^bMS /A^rMS ratios by extracting radial bins 
using simple circles, and we found no significant differences be- 
tween the two radial binning methods. 

Estimating the eiTors of our points required special attention. 
First we took the Poisson errors of the numbers of stars, and used 
them to generate Poisson errors for the values of A^bMS /A^rMS ■ 
These, however, are only a lower bound for the true error, which 
has additional contributions that are impossible to estimate di- 
rectly; they come from blends, binaries, etc. To estimate the true 
eiTors empirically, for each value of A^bMs/Mivis we subdivided 
the sample of stars that had been used. In the inner ACS/WFC 
mosaic the subsamples were the quadrants shown in Fig.[T] while 
for each of the outer fields, where we do not have symmetric 
azimuthal coverage, we divided the sample into magnitude in- 
tervals, four for each FORSl @VLT field and three each for the 
WFPC2 field and the outer ACS/WFC field. 

We treated each set of subsamples as follows: Within each 
subsample we performed a dual-Gaussian fit, and derived from 
it the value of A^bMs/A^rMS- We weighted each subsample accord- 
ing to the number of stars in it, and took a weighted mean of 
the four (or three) values of A^bMs/M-MS, to verify that this mean 
was equal, within acceptable round-off errors, to the value that 
we had found for the whole sample. (It was, within a per cent 
or two in nearly every case.) Finally we derived an error for the 



sample, from the residuals of the individual A^bMS /Mms values 
from their mean, using the same weights as we had used for 
the mean. These are the error bars that are shown in Figure [10] 
These errors are indeed larger than the Poisson errors, but only 
by about 10%. We must note, however, that in addition to the 
random error represented by the error bars, it is likely that there 
is still some systematic eiTor in our values of A^bMs/A^rMS, due 
to the effects of blends and binaries. On the one hand, blends 
have the same photometric effect as true binaries; they tend to 
move bMS stars into the rMS region, while many of the rMS 
stars that are similarly affected are eliminated by our red cut-off. 
This effect tends to reduce our observed value of A^bMs/Mivis- It 
is less easy to predict, however, how such effects increase toward 
the cluster center. Blends, on the one hand, increase because of 
the greater crowding. Binaries, on the other hand, increase be- 
cause their greater mass gives them a greater central concentra- 
tion. To repeat, the result has been that our values of A^bMs/A^rMS 
are somewhat depressed toward the cluster center, so that the 
gradient of A^bMS /Mms that we report is probably a Httle lower 
than the real one. 

Fig. [To] clearly shows a strong radial trend in the ratio of 
bMS to rMS stars, with the bMS stars more centrally concen- 
trated than the rMS stars. The most metal-rich population, MS- 
a, is too sparse, and also too hopelessly mixed with the red edge 
of the rMS, to allow any reliable measurement of its radial dis- 
tribution, but in the next section we will examine the distribu- 
tion of its progeny, RGB-a. Table [2] summarizes our results. The 
first column identifies the data set. Columns 2-4 give, for the 
inner ACS/WFC 3x3 mosaic, the minimum, median, and max- 
imum radius of the central circle or the annulus, while for the 
other fields these columns give the inner, median, and maxi- 
mum radius that the field covers. The sigmas of the Gaussians 
that best fit the bMS and rMS color distributions, with their 
uncertainties, are in Columns 5-8. Columns 9 and 10 give the 
A^bMS /A^rMS ratio and its eiTor. Column 1 1 gives the difference 



12 



Bellini, A. et al.: Radial distribution of the multiple stellar populations in oj Centauri 



(in straightened color) between the peaks of the Gaussians that 
best fit the bMS and rMS. The last column gives the color base- 
line of each data set. By Acolor we mean a color difference or 
width, in the straightened CMD [either (B, B - Rc) or (otf606w, 
wipeoew - 'Wf814w), whichever applies]. 

Our results are qualitatively consistent with those of Sollima 
et al. (120071) . within the common region of radial coverage. We 
confirm the flat radial distribution of A^bMs/Mivis outside ~8-10 
arcmin, and a clear increase of A^bMs/Mivis toward the cluster 
center. For the first time, and as a complement to the Sollima 
et al. (I2007I I investigation, our ACS/WFC 3x3 mosaic data set 
has enabled us to study the distribution of ll) Cen MS stars in the 
innermost region of the cluster Inside of ~1.5 rc (i.e., inward 
of ~2'), the A^bMS /A^rMS ratio is almost flat and close to unity, 
with a slight overabundance of bMS stars. At larger distances 
from the cluster center, the A^bMS /A^rMS ratio starts decreasing. 
Between ~3' and ~8' (the latter corresponding to ~2 half-mass 
radii) the ratio rapidly decreases to ~0.4, and remains constant 
in the cluster envelope. Better azimuthal and radial coverage of 
the region where the maximum gradient is observed would be 
of great value. In the radial interval between 1 and 2 half-mass 
radii, we can use only the corners of the ACS/WFC 3x3 mo- 
saic, and the FORSl photometry, which inside of 10' is seri- 
ously affected by crowding and saturated stars. In any case, the 
star counts and even visual inspection of the histograms in Fig.|9] 
leave no doubt about the overall gradient. 

Note that in the two innermost bins there are more bMS than 
rMS stars, even though the heights of the two peaks would sug- 
gest the opposite. The apparent contradiction disappears, how- 
ever, when we note the much greater width of the bMS Gaussian, 
which more than makes up for the difference in heights. This 
seems to be consistent with a greater spread in chemical com- 
position for metal-intermediate than for metal-poor stars, as 
first seen by Norris et al. ( 11997 1. Our approach, using a dual- 
Gaussian fit, has been optimized to estimate the value of the 
number ratio of bMS to rMS stars, avoiding as much as possible 
any contamination by blends, binaries, and MS-a stars. 

We must also address the fact that the A^bMs/A^rMS values 
found by Sollima et al. ( 120071 ) are consistently lower than our 
values. The difference is largely due to their use, on the red side, 
of a wide color range (see their Fig. 5) that includes nearly all 
of the contamination by blends, binaries, and MS-a stars that 
our method has so studiously avoided. This makes their num- 
bers of rMS stars much too high — easily enough to account for 
their finding a value of ~0. 16 in the cluster envelope, rather than 
our ~0.4, which is certainly much closer to the truth. Note also 
that we have concentrated exclusively on the ratio of numbers of 
bMS and rMS stars, making no attempt to derive absolute num- 
bers for each component. We felt that absolute numbers would 
be subject to different incompleteness corrections in our differ- 
ent data sets, whereas the incompleteness in each data set should 
be the same for each component and should therefore not affect 
their ratio. 

Finally, the robustness of our method is shown by the close 
agreement of our — proper-motion selected — outer ACS field 
(magenta open circle in Fig.fTOli. which has almost no crowding 
problems, with the outer ground-based FORS 1 fields (last two 
red squares in Fig. [TOl ), which are certainly affected somewhat 
by crowding. 

3.4. Artificial star tests 

Even with the technique that we have used to exclude the ef- 
fects of photometric blends and binaries, which lie above and 



to the red of any MS, there is a concern that some bMS stars 
would be shifted into the rMS region (and some rMS stars lost 
on the red side of the MS), and that these shifts would distort 
the NhMs/NrMS ratio. As a check against this possibility we have 
made two tests using artificial stars (AS). In each test we intro- 
duced the same AS into both the F435W and the F625W images, 
as follows. 

For each test, we first created 45 000 artificial stars, with ran- 
dom F435W instrumental magnitudes between -11.1 and -9.9 
(corresponding to 20.9 < B < 21.1), and random positions. We 
then took each of these 45 000 AS, assigned it a color that placed 
it on the bMS, and inserted it in the F625W images, at the same 
position but with the F625W magnitude that corresponds to this 
color We then repeated this procedure for 45 000 new AS, but 
this time we gave them colors that put them on the rMS. (What 
we mean by "on" [bMS or rMS] differs between the two tests; 
see below for an explanation of the difference.) Each artificial 
star in turn was added, measured, and then removed, so as not 
to interfere with the other AS that were to be added after it; this 
procedure is that of Anderson et al. (2008), where it is explained 
in detail. 

In order to test the effects of crowding, each of the two 
tests used two fields from the central 3x3 mosaic: the central 
field where crowding effects are maximal, and one of the corner 
fields, about 5 arcmin (3.6rc) southeast of the center (see Fig.[T] 
for a map of the 3x3 mosaic of fields). 

The first of the tests (TESTl) was aimed at checking the pho- 
tometric errors in the colors. To do this, we chose the color of 
each AS so as to put it exactly on the ridge line of the bMS or 
rMS; the color spread of the recovered AS would then serve as a 
lower-limit estimate of our photometric error. 

The aim of TEST2 is to verify our ability to insert AS with 
A^bMs/A^rMS=l and then recover that value, when the two MSs 
have intrinsic dispersions in color. To do this, we first derived 
the intrinsic spreads of the two sequences by taking from the fifth 
and seventh columns of Table |2] the simple unweighted mean of 
the entries in lines 1 and 2 for the central field, and in lines 4 and 
5 for the corner one. (The more fastidious procedure, weighting 
the entry in each of the two lines according to the number of stars 
contributed by that annulus, would have been quite laborious and 
would have made no significant change in the results.) These are 
the observed total color spreads (intrinsic spread plus measur- 
ing error) of the bMS and rMS, respectively, in the two fields 
that we are using here. From these total spreads we quadrati- 
cally subtracted the corresponding measuring-error spreads that 
we had found in TESTl, so as to get estimates of the intrinsic 
color spreads of the two sequences. We created AS in the same 
manner as in TESTl, but this time instead of placing the AS on 
the center lines of bMS and rMS, we adjusted the F625W mag- 
nitudes so as to give the AS a Gaussian spread in color around 
each sequence, using the intrinsic sigmas that we had just found. 
After the measuring process, these AS should duplicate the ob- 
served total spreads, and can be used to estimate the amount of 
contamination between the two main sequences. To repeat, each 
test was performed both on both the central and the corner field. 

The results of these AS tests are summarized in brief numeri- 
cal form in Tableland in graphical form in Figures[TT|and[T2] In 
each figure the left and right halves refer to the central and corner 
fields, respectively, while each half figure is divided into three 
panels that show, from left to right, the CMD, the straightened 
CMD, and the decomposition of the number densities of the lat- 
ter into best-fitting Gaussians. Each panel showing the Gaussian 
fits is subdivided into five magnitude intervals, (very similarly 
to what is done in Fig. 8). Cols. 2-4 of Table [3] give, for each 
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Table 3. Results of the two artificial-star tests. For each of the two fields (first column), we give in Cols. 2^ the values of A^bMs/Misis 
for the AS that were inserted, and the color dispersions that were given to the AS that were put on the bMS and rMS, respectively. In 
Cols. 5-8 are, respectively, the A^bMS /A^rMS of the AS that were recovered, with error, followed by the sigmas of the two Gaussians 
that were fitted to them. See text for details. 



Field (A'bMs/-'VrMs)ins. (0"bMs)ins. (0"rMs)ms. (^bMS /-'ViMS )rcc. 0"(^bMS /A'rMs)rec. (0"bMs)rcc. (0"rMs)rcc. 



TESTl 



central 
comer 



1.000 
1.000 



0.000 
0.000 



0.000 
0.000 



0.946 
0.992 



0.013 
0.013 



0.013 
0.009 



0.013 
0.009 



TEST2 



central 
comer 



1.000 
1.000 



0.021 
0.016 



0.016 
0.013 



1.008 
0.996 



0.079 
0.027 



0.026 
0.019 



0.021 
0.016 



TESTl (central field) TESTl (corner field) 




Fig. 11. TESTl artificial star analysis for the central ACS/WFC 3x3 mosaic field (left panels), and for the corner field (right panels). 
For each panel, we show the CMD with the recovered stars (in blue for the bMS stars and in red for the rMS stars), for five magnitude 
intervals. The straightened MSs are plotted in the middle, while on the right we show the color histograms, with the dual-Gaussian 
fits. The vertical lines in blue, red, and black mark, respectively, the centers of the two Gaussians and the red cut-off. See text for 
details. 



TESTS (central field) TEST3 (corner field) 




Fig. 12. Same as Fig.[TTl but now for TEST2. See text for details. 



field and AS test, the A^bMs/Mivis ratio of the inserted AS and the of the five magnitude bins and its error, as explained in detail for 
dispersions of the MSs. The recovered values (weighted mean real stars in Sect. 3.3) are shown in Cols. 5-8. 
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From the results of TESTl we conclude that in each field 
the color spread introduced by measuring error is the same for 
bMS stars as for the rMS, and that it is about 40% higher in the 
central field than in the less-crowded corner field. TESTl has 
served two purposes: (1) It gave us a clear, effective measure 
of the effect of crowding on the color spread. (2) It evaluated 
the color spreads due to measuring error alone, which we used 
in setting up TEST2. (Its results for A^bMs/M-MS are given, pro 
forma, but they have no real significance, since the color spreads 
used in TESTl are so naiTow that our color bin-width does not 
sample them adequately.) It is TEST2 which directly tests our 
previous conclusions about the size of A^bMs/M-MS- We conclude 
from it that the AS tests recover our input values of A^bMs/A^rMS, 
within the uncertainties of the measurement. 

In this section we have demonstrated, on two extreme fields 
of the ACS inner mosaic, that our dual-Gaussian fitting method 
is fully effective in overcoming the effects of crowding on the 
distribution of colors, and that it reliably estimates the relative 
star numbers in the two sequences. (Note that we use this same 
method for all of our other data sets too). As we noted at the end 
of Sect. 3.3, the excellent agreement between the results from 
our completely uncrowded outer ACS field and those from our 
outer FORS 1 fields establishes the validity of the latter, without 
recourse to any additional AS texts for them. 

4. Radial gradients in the RGB subpopulations 

It has been known since the end of the 60s that the RGB of 
oj Cen is broader than would be expected from photometric er- 
rors (Woolley & Dickens 1967), but it was only in 1999 that Lee 
et al. ( 1999) clearly detected at least two distinct RGBs. Later on, 
Pancino et al. ( I2000I I demonstrated that there is a correlation be- 
tween the photometric peaks across the RGB and three peaks in 
the metallicity distribution. On this basis, they defined the three 
RGB groups: RGB-MP, RGB-Mint, and RGB-a, characterized 
by an increasing metallicity. In this section we will present a de- 
tailed study of the radial distributions of these components. 

4.1. Defining fhe RGB-MP, RGB-Mint, and RGB-a 
subsamples 

Unfortunately the WFPC2, FORSl, and outer ACS/WFC data 
sets we used to analyze the main-sequence population in the pre- 
vious section are saturated even at the MS turn-off level, and are 
therefore unusable for the study of the RGB radial distributions. 

Our WFI@2.2m photometric and proper- motion catalog 
(Bellini et al. 120091 ), however, is an excellent data base for this 
study, particularly in view of the fact that we can safely remove 
field objects in the foreground and background, thanks to our 
accurate proper motions. This proper-motion cleaning is of fun- 
damental importance in the outer envelope of the cluster, where 
there can be more field stars than cluster giants. In the cen- 
tral regions of the cluster, the WFI@2.2m data are less accu- 
rate due to the poorer photometry caused by the crowded con- 
ditions, so there we take advantage of our high-resolution inner 
ACS/WFC 3x3 mosaic, which included short exposures to mea- 
sure the bright stars. Below we describe how we extracted the 
(L> Cen RGB subsamples from these two data sets. 

Because of the complex distribution of the stars along the 
RGB we were forced to use bounding boxes to select the differ- 
ent RGBs. This selection procedure is less accurate than what we 
were able to do for the bMS and the rMS; nevertheless it is still 
accurate enough to study the general trend of the radial distribu- 
tion of the relative numbers of RGB-MP, RGB-Mint, and RGB-a 



stars. The Poisson error from the smaller number of RGB stars 
makes the more precise procedure less critical. 

For the ACS data, we defined bounding boxes for the RGB 
subpopulations of oj Cen in the CMD obtained from the data 
set of the ACS/WFC 3x3 mosaic, for which the large-number 
statistics make the separation among the different RGBs easier 
to see. We extracted three RGB subpopulations, in a way very 
similar to that used by Ferraro et al. (.2002) . [Note that other 
authors (e.g. Rey et al. 120041 SoUima et al. |2005al Johnson et 
al. 120091 ) have defined four or even five RGB subpopulations]. 
The left panel of Fig. [T3] shows the three RGB bounding-box 
regions drawn in the CMD from the ACS/WFC 3x3 mosaic, to 
identify the three subgroups RGB-MP, RGB-Mint, and RGB- 
a. Our RGB selections are limited to magnitudes brighter than 
B = 17.9, and contain 5184 RGB-MP stars, 4379 RGB-Mint 
stars, and 383 RGB-a stars. 

In extracting the RGB subpopulations from our WFI@2.2m 
data set we chose to define the subpopulations in the B, B - V 
CMD. Even though we cannot adopt exactly the same selection 
boxes in the B, B - Rc CMD as for the ACS/WFC 3x3 mosaic. 
This choice might appear awkward, not only because the color 
baseline B - V is shorter than the B - Rc baseline, but also be- 
cause the WFI Rc filter is very similar to the ACS/WFC F625W 
filter. There are other good reasons for adopting the B - V color 
baseline, however. The most important one is that the WFI pho- 
tometry obtained with the V filter has ten times as much integra- 
tion time, and more dithered images than those available for the 
Rc filter. Therefore our V photometry is considerably more pre- 
cise, and more accurate, than our Rc magnitudes. Moreover, our 
empirical sky-concentration correction (very important for such 
studies) is better defined in V than in Rc (see Bellini et al. 120091 ). 

In this WFI@2.2m B vs. B - V CMD, we tried to define the 
bounding boxes in a way that was as consistent as possible with 
what we did for the data set from ACS/WFC 3x3 mosaic. We 
cross-identified the stars that are in common between the sample 
that we had selected from the RGB CMD of the ACS/WFC 3x3 
mosaic, on the one hand, and the WFI@2.2m B - V data set on 
the other hand, and we carefully drew by hand, in the (B, B -V) 
CMD, bounding boxes that would include the same stars as in 
the sample from the ACS/WFC 3x3 mosaic. 

In addition, we selected from the WFI@2.2m data set the 
stars that were measured best (both photometrically and astro- 
metrically), and were most likely to be members of a> Cen. To 
make the selection we used the eiTor quantities in columns 7, 9, 
13, and 15 of Table 6 of Bellini et al. (2009) . These are the er- 
rors of the two components of proper motion and of the B and 
V magnitudes. Our selection consisted of choosing, at the bright 
end of the RGB, stars whose proper-motion error has a magni- 
tude less than 1.8 mas yr"', and whose photometric eiTor is less 
than 0.02 mag in each band,; we also required that the proper 
motion of a star differs from the mean motion of cluster stars 
by no more than 2.1 masyr At the faint end of the RGB we 
allowed these three tolerances to rise to: 2.1 masyr"', 0.03 mag 
and 3.8 masyr"', respectively. This high-quality data set com- 
prised 4993 RGB-MP stars, 3057 RGB-Mint, and 292 RGB-a 
stars. 

The right-hand panel of Fig. [T3]shows the WFI@2.2m RGB 
subpopulations that were selected in this way. We note that 
whereas the RGB-a sample is well separated from the other two 
RGB components, the RGB-MP and RGB-Mint components are 
separated only by an arbitrary dividing line, so that small dif- 
ferences in defining the bounding boxes might result in some 
cross-contamination in those two samples. 
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B-R^ 



B-V 



Fig. 13. CMDs of the u Cen RGB from ACS/WFC 3x3 mosaic 
{B vs. B-Rc, left panel) and from WFI@2.2m data {B vs. B-V, 
right panel). The RGB subpopulations selected are also plotted 
with different colors. See text for details. 




Fig. 14. (a): Radial distribution of the ratio RGB-a/RGB-MInt 
for the WFI@2.2ot data set (red triangles) and for the ACS/WFC 
3x3 data set (blue circles). Vertical dashed lines mark the core 
radius and the half-mass radius, respectively, (b): Radial distri- 
bution of the ratio RGB-MInt/RGB-MR (c): Radial distribution 
of the ratio RGB-a/RGB-MR See text for details. 



4.2. Relative radial distributions of RGB stars 

We divided our WFI@2.2m data set into ten radial bins, each 
containing approximately the same number of RGB-Mint stars, 
and the ACS/WFC 3x3 data set into five radial bins, again with 
the same equal-number criterion. For each of these bins we 
counted the number of RGB stars in each subpopulation. 

In Fig.[T4]we show the derived radial gradients. As it has not 
been possible to perform the same error analysis as was done for 
the MS stars (because of the much smaller number of stars), the 
error bars in Fig.fT4li'epresent only Poisson eiTors, and should be 
considered a lower limit to the real errors. In panel (a) we show 
the radial distribution of the ratio RGB-a/RGB-MInt. Blue full 
circles refer to the ACS/WFC 3x3 data set, and red triangles to 
the WFI@2.2m data. Vertical dashed lines mark the core radius 
rc and the half-mass radius r/,. We found that, within the eiTors, 
the RGB-a and the RGB-Mint stars share the same radial dis- 
tribution, since their ratio is constant over the entire radial range 
covered by our two data sets. In panel (b), we plot the ratio RGB- 
Mlnt/RGB-MP for the two data sets. The RGB-Mint stars are 
more centrally concentrated than the RGB-MP stars, with a flat- 
ter trend within ~1 r/,, a rapid decline out to ~8'-10', and again 
a flat relative distribution outside. There is a hint, also, that the 
RGB -Mint/RGB -MP ratio could be nearly constant within the 
half-mass radius. We find that the general radial trend of the 
RGB -Mint/RGB -MP star-count ratio is consistent with that of 
M)MS /^rMS • This result provides additional evidence (in agree- 
ment with the metallicity measurements by Piotto et al. 120051) 
that the bMS and the RGB-Mint population must be part of the 
same group of stars, with the same metal content and the same 
radial distribution within the cluster. Panel (c) shows that the 
ratio RGB-a/RGB-MP resembles, within the errors, the RGB- 
Mlnt/RGB-MP trend. We were unable to examine this trend for 
the MS part of the RGB-a population, since the MS-a sequence 
cannot be followed below B ~20. 

Our analysis confirms the results by Norris et al. ( 1 19971 1. 
Hilker & Richtler (TOOO), Pancino et al. (l2000b . and Rey et 
al. (2004), and Johnson et al. (2009) , who found that the most 
metal-poor RGB stars are less concentrated than the RGB-Mint 
ones. Moreover, we can also confirm that the RGB-a and the 
RGB-Mint share the same radial distribution within co Cen, as 
found by NoiTis et al. ( 1997), Pancino et al. ( 120001 ). and Pancino 
et al. (2003) for RGB-a only. 

It is important to note that because we were able to use 
proper motions to construct a pure cluster sample, our results 
are not affected by field-star contamination, which would tend 
to enhance the RGB-a star counts in the cluster outskirts with 
respect to the more populous RGB-MP sample (which also cov- 
ers a smaller region in the CMD). Field-star contamination is 
likely the reason that Hilker & Richtler ( 120001 ) and Castellani et 
al. (2007 ) found the RGB-a/RGB-MP ratio to increase with dis- 
tance from the cluster center — the opposite trend from what 
is seen here. Moreover, it is interesting to note that the dif- 
ferent RGB-Mint subgroups (as highlighted, e.g., by Sollima 
et al. |2005al ) might well have a different radial behavior, but 
necessarily — since we cannot distinguish them in the CMD — we 
have to treat them together and study only their average gradient. 



5. Discussion 

In this paper we have analyzed the radial distribution of the 
different MS and RGB components in the globular cluster a> 
Centauri. We used high-resolution ACS/WFC images to study 
the inner regions of the cluster, and ACS/WFC, WFPC2 and 
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FORSl@VLT images, as well as WFI@2.2m images, for the 
cluster envelope. We found that there are slightly more bMS 
stars than rMS stars in the inner 2 core radii. At larger distances 
from the cluster center, out to ~8 arcmin, the relative number 
of A^bMS /Mms stars drops sharply, and then remains constant at 
M)Ms/A^rMS~0.4, out to half the tidal radius of the cluster 

Our most precise photometry comes from the outer ACS 
field at 17' (12 r^), where we find that the color dispersion (cr) 
of the bMS is about 50% larger than that of the rMS. The other 
observations are consistent with this, though they are unable to 
measure cr so precisely, on account of crowding (in the inner 
ACS field) and other errors (in the ground-based fields). 

The R GB-Mi nt population (associated with the bMS by 
Piotto et al. l2005l l and the RGB-MP sample (which includes the 
progeny of the rMS) follow a trend similar to that of A^bMS /A^rMS ■ 
The most metal-rich component of the RGB, RGB-a, also fol- 
lows the same distribution as the RGB-Mint component. 

On the hypothesis that the bMS, the presumably helium-rich 
population, is a second generation of stars formed by the low- 
velocity material ejected by a primordial population (which we 
assume to be the more metal-poor rMS population), the bMS 
must have formed from matter that collected in the cluster cen- 
ter via some kind of cooling flow. This is in qualitative agree- 
ment with the recent models by Bekki & Norris (120061) and 
D'Ercole et al. (2008). The very long relaxation time (half- 
mass relaxation time longer than 10 Gyr, according to the Harris 
119961 compilation) has preserved some information about the 
original kinematic and spatial distribution of the material from 
which the younger component took form. Interestingly enough, 
the third, most-metal-rich population is also more concentrated 
than the most metal-poor component, and has a radial distribu- 
tion that is rather similar to that of the intermediate-metallicity 
sample. It is also noteworthy that the bMS component has a 
broader color distribution than the rMS one. This fact may re- 
flect, at least in part, the large dispersion in iron abundance of 
the intermediate-metallicity component (e.g. Norris & Da Costa 
I1995I I. Alternatively, this bMS spread could be an indication of 
the dispersion of other chemical elements, including He. Only a 
detailed analysis of the metal content of the two MSs can solve 
this issue, but for this we might need to wait for the next gen- 
eration of 30+ meter telescopes, on account of the faintness of 
these stars. 
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